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Abstract. We present the results of numerical calculations of the channelling and 
undulator radiation generated by an ultra-relativistic positron channelling along a 
crystal plane, which is periodically bent. The bending might be due either to the 
propagation of a transverse acoustic wave through the crystal, or due to the static 
strain as it occurs in superlattices. The periodically bent crystal serves as an undulator. 
We investigate the dependence of the intensities of both the ordinary channelling 
and the undulator radiations on the parameters of the periodically bent channel with 
simultaneous account for the dcchannclling effect of the positrons. We demonstrate 
that there is a range of parameters in which the undulator radiation dominates over 
the channelling one and the characteristic frequencies of both types of radiation are 
well separated. This result is important, because the undulator radiation can be used 
to create a tunable source of X-ray and 7-radiation. 



In this Letter we present for the first time the results of calculations of the total 
spectrum of emitted photons accompanying the channelling process of ultra-relativistic 
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charged particles through a crystal which is periodically bent. The present consideration 
is a further step in investigating new phenomenon which was described recently in 

|2| and was called Acoustically Induced Radiation (AIR). It was noted that the 
periodic pattern of crystal bendings (which can be achieved either through propagation 
of a transverse acoustic wave, AW, or by using static periodically strained crystalline 
structures p), [|) gives rise to a new mechanism of electromagnetic emission of the 
undulator type, in addition to a well-known ordinary channelling radiation ||. 

Without any loss of generality further in this Letter we consider the case of a 
dynamic periodic bending of a crystal by means of transverse AW, as it is illustrated 
in figure [p. Under the action of the transverse AW propagating along the ^-direction, 
which defines the center line of initially straight channel (not plotted in the figure) 
the channel becomes periodically bent. Provided certain conditions are fulfilled || the 
beam of particles, which enters the crystal at a small incident angle with respect to the 
curved crystallographic plane, will penetrate through the crystal following the bendings 
of its channel. This results in transverse oscillations of the beam particles (additional to 
the oscillations inside the channel due to the action of the interplanar force). These 
oscillations become an effective source of spontaneous radiation of undulator type 
due to the constructive interference of the photons emitted from similar parts of the 
trajectory. It was demonstrated that the system "ultra-relativistic charged particle 
+ periodically bent crystal" serves as a new type of undulator, and, consequently, as a 
new source of undulator radiation of high intensity, monochromaticity and of a particular 
pattern of the angular-frequency distribution. 



spontaneous undulator 
emission 
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Figure 1. Schematic representation of spontaneous and stimulated AIR in a channel 
bent by the transverse AW. The y- and z-scales are incompatible. 



As it was pointed out in |2j this scheme leads, in addition to the spontaneous 
radiation, to a possibility to generate stimulated emission, similar to the one known 
for a free electron laser |5| in which the periodicity of a trajectory of an ultra-relativistic 
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projectile is achieved by applying spatially periodic magnetic field. In connection with 
the stimulated AIR it was noted in that to achieve noticeable degree of amplification 
one has to operate with a positron bunch of a high volume density, which, nevertheless, 
turned out to be achievable in modern accelerators as discussed in ||. 

The main subject of our previous studies |2j was the characteristics of the 
undulator AIR itself. Due to this reason we primarily investigate the case of a high- 
amplitude AW, a^> d with d standing for the interplanar spacing, propagating through 
the crystal. Less attention, except for some qualitative estimates, has been paid to the 
detailed investigation of the mutual influence of two types of radiation, the ordinary 
channelling radiation and the AIR, in forming the total spectrum of the radiation 
emitted by a channeled positron. This topic is addressed in the present Letter where 
we report our first quantitative results on numerical calculation of the spectrum of 
the emitted photons with both mechanisms taken into account simultaneously. We 
demonstrate that there are ranges of (i) the parameters of AW, which are the amplitude, 
a, the wavelength A, and the sound velocity V, (ii) the energies of projectile positron, e, 
(iii) the crystal parameters, which include the length of a crystal, the constituent atoms 
and the types of crystallographic planes, inside which 

• the characteristic frequencies of the AIR and the ordinary channelling radiation are 
well separated, 

• the intensity of AIR is essentially higher than that of the ordinary channelling 
radiation, 

• the radiative spectrum is stable towards the total losses of the particle (in the case 
of a positron it is primarily the radiative ones) 

These items, except the last one which was considered in J7J, are discussed below in 
the Letter. 

To conclude the introductory part we mention that in the past the problem of 
evaluating the total spectrum of radiation formed in a bent crystal was considered in 
several publications || [| in the case of a projectile channelling in a crystal bent with a 
constant curvature radius. In our Letter for the first we investigate the problem for the 
periodically bent channel. In full the results of our research will be published elsewhere 
|10| . Below we present the essential points. 

An adequate approach to the problem of the radiation emission by an ultra- 
relativistic particle moving in an external field was developed by Baier and Katkov 
in the late 1960s [1 1| and was called "operator quasi-classical method" by the authors. 
The details of that formalism can be found in [ 12| , |T3|j . The advantage of this method 
is that it allows to use the classical trajectory for the particle in an external field and, 
simultaneously, it takes into account the effect of the radiative recoil. 

For particles with spin s — 1/2 the energy radiated into a given direction n summed 
over the polarizations of the photon and the projectile is given by (the CGS system is 
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adopted throughout the paper) 

where a ~ 1/137 is the fine structure constant, <^(ti,t 2 ) = ti — h — n- (r(ii) — r(t 2 ))/c, 
and /(ti,t 2 ) = {[1 + (1 +uf] (v(t!)v(t 2 )/c 2 - 1) +u 2 7~ 2 } /2, where c is the velocity of 
light, 7 = e/m is the relativistic factor and u — Two/{e — TuS). 

The main advantage of the the quasi-classical approach is that to calculate 
the angular-spectral distribution of the radiation one only needs to know the time 
dependencies of classical radius-vector r(i) and the velocity v(i) of projectile. In 
connection with a positron channelling through a periodically bent crystal the above 
relations impose, if applied directly, some restrictions on the projectile energy, on the 
parameters of the crystal (channel) and the AW. 

The channelling process in a bent crystal takes place if the centrifugal force in the 



channel is less than the maximal force due to the interplanar field |14| . For a periodically 
bent crystal the maximal centrifugal force is equal to rri'yv 2 /R min , with v ~ c and -R m i n 
being a minimum curvature radius of the bent channel. Hence, the following condition 
must be fulfilled |L], § 

m 7 c 2 /i? min < U' max . (2) 

where U' max stands for the maximum gradient of the interplanar field. For an acoustically 
bent channel -R m i n = (A/27r) 2 /a, therefore, the inequality (0) bounds together the 
characteristics of the crystal, £7 max , the AW amplitude and the wavelength, and the 
relativistic factor 7. A comprehensive study of the allowed ranges of all parameters 
involved in @ was carried out in [|l], [|. 

Second requirement which has to be fulfilled to make eq. ([!]) directly applicable to 
calculating the spectrum accompanying the channelling process concerns the upper limit 
of integration r, which is related to the crystal length L through L = ct, and which, 
for given value of the AW wavelength, defines the number of the undulator periods as 
N = L/X. The AIR acquires specific features of the undulator-type radiation (such as 
the monochromaticity and particular pattern of the angular-frequency distribution, see 
e.g. |l!J) provided N ^> 1. The length of a crystal is subject to a physical condition that 
the positron bunch stays inside the channel when penetrating through the crystal on the 
scale of L. In reality, the parasitic effect, the dechannelling leads to a decrease in the 
volume density of the channeled particles n(z) with penetration distance z, and roughly 



satisfies the exponential decay law for both straight and bent channels (see |15|, 
n(z) = no exp (—z/Ld), where no is the volume density t the entrance, and La is the 
dechannelling length, which for given crystal and channel depends on a positron energy 
and on the curvature radius. For a periodically bent crystal Ld(j, R) can be estimated 
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as g 0, 



L rf(7 ,i?) = (l-^ mm ) 2 ^^^7 (3) 



where r c i is the classical radius of an electron, otf = 0.8853Z _1 / 3 ao (ao is the Bohr 
radius) and 7 = 16Z 9 eV are, respectively, the Thomas-Fermi radius and ionization 
potential of the crystal atoms, Z is the atomic number, d is the interplanar distance in 
the lattice. The quantity R c = £/U^ nax is the critical (minimal) radius consistent with the 
channelling condition in a bent crystal (@). The quantity L c = In (V^ymc 2 //) - 23/24 
is the Coulomb logarithm characterizing the ionization losses of an ultra-relativistic 
positron in amorphous media with account for the density effect (see e.g. [[[7], pp. 

The AIR will have a pronounced pattern of undulator-type radiation provided the 
number of the periods is large on the scale of Combined with (|3]) this condition 
leads to another restriction on the parameters involved 

where A can be related to the AW frequency v through A = V/v with V standing for 
the sound velocity. 

Figure |2| illustrates the restrictions which are imposed on the values of a and v 
by inequalities (0) and (f|) in the case of e — 0.5 GeV planar channelling in Si along 
(110) crystallographic planes. The diagonal straight lines correspond to various values 
(as indicated) of the parameter C = £ / (RmmU^x) — 1 consistent with the channelling 
condition (^). The curved lines correspond to various values (as indicated) of the number 
of undulator periods N related to the dechannelling length L d through eq. (f|). The 
horizontal lines mark the values of the AW amplitude equal to d (with d = 1.92 A being 
the (110) interplanar distance in Si) and to 10 d. The vertical line marks the value 
v = 200 MHz of the AW frequency for which the spectra ([I]) were calculated. We used 
the value V = 4.67 x 10 5 cm/s for the velocity of sound in Si (this value was obtained 
by using the data from [0). Thus, the AW wave-length used in our calculations equals 
to A = 2.33 x 10~ 3 cm. 

The calculated spectra of the radiation emitted in the forward direction (with respect 
to the z-axis, see figure P for photon energies from 45 keV to 1.5 MeV are presented 
in figures [| The details of the analytical evaluation of the right-hand side of (HD and 
its numerical implementation as well as more extended results of the calculation will 
be published soon ||10|| . Here we briefly sketch the numeric procedure and present a 



discussion of the exhibited results. The AW frequency, the number of undulator periods 
and crystal length were fixed at v = 200 MHz, iV = 15 and L = NX = 3.5 x 10~ 2 cm. 
The ratio a/d of the AW amplitude to the interplanar spacing was varied within the 
interval a/d = ... 10. The case a/d = corresponds to the straight channel. 
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100 

AW frequency, v (MHz) 

Figure 2. The ranges of the allowed parameters of the AW, a versus v , consistent 
with the conditions (||) and for e = 0.5 GeV (7 s» 10 3 ) positron channeling in Si 
along (110) crystallographic planes. See explanations in the text for the meaning of 
various lines. 



To evaluate the spectral distribution ([[]) the following procedure was adopted. 

Firstly, the spectrum was calculated for individual trajectories. These were obtained 
by solving the relativistic equations of motion with both the interplanar and the 
centrifugal potentials taken into account. We used the continuum approximation [^(| 
to describe the projectile positron - lattice atoms interaction. Within these scheme 
we considered two frequently used |]15| analytic forms for the interplanar potential, 
the harmonic and the Moliere potentials calculated at the temperature T = 150 K to 
account for the thermal vibrations of the lattice atoms. For each a/d value by changing 
the initial values of the entrance coordinate y^ and the initial velocities along the y-axis 
(see figure [[]) we left only those trajectories which corresponded to the case of stable 
channelling through the whole crystal length L. We call a trajectory as a "stable" one 
if moving along it the particle does not approach crystalline planes at a distance less 
than the Thomas- Fermi radius <2tf (^tf = 0.194 A for a Si atom). This allowed us 
to totally disregard, at least on the scale L ~ La, the random scattering of a projectile 
by lattice electrons (see e.g. |15|, |l6|j). Thus, for each a/d value, we defined the ranges 
of the initial coordinates G [— d/2 + a TF ,c?/2 — a TF ] and the velocities vf \ and, 
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Figure 3. Spectral distribution (JsJ) of the total radiation emitted in the forward 
direction (n||z) for e — 0.5 GeV (7 » 10 3 ) positron channeling in Si along (110) 
crystallographic planes calculated at different a/d ratios as indicated. The AW 
frequency is fixed at v = 200 MHz, the crystal length is L = 3.5 x 10~ 2 cm. 



correspondingly, the initial phase volume &°'(a/d) = §p^dy^ (where stands 
for initial transverse momentum) for which the corresponding classical trajectories are 
stable when channelling through the whole crystal length L. Then, discretizing the found 
initial phase volume <$>(°\a/d) by choosing N y (o) x iV(o) points (?/ (0) ,^ 0) ) e $(°)(a/d), 
the individual spectra dE(y(°\p^) were calculated for each pair of the initial coordinate 
and velocity. Finally, for each a/d value we evaluated the averaged spectrum defined as 
follows: 

Here, the integration is carried out over the phase volume &(°'(a/d), and the integral 
is scaled by the phase volume &°>(a/d = 0) of stable trajectories in straight channel. 
The ratio (a / d) / (a / d = 0) describes the number of particles channelled through 
the bent crystal relative to the number of particles channelled through the straight one. 
Hence, it is convenient to use the quantity (dE) to compare the spectra produced by 
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effectively different number of projectiles as it happens for different a/d values. 

Figures [| correspond to the spectra (Q) obtained as outlined above. The results 
presented were calculated by using the Moliere approximation for interplanar potential. 

The first graph in figure ^| corresponds to the case of zero amplitude AW (the ratio 
a/d = 0) and, hence, presents the spectral dependence of the ordinary channelling 
radiation only. The asymmetric shape of the calculated ordinary channelling radiation 
peak, which is due to the strong anharmonic character of the Moliere potential, bears 
close resemblance with the experimentally measured spectra The spectrum starts 
at hu ~ 960 keV, reaches its maximum value at 1190 keV, and steeply cuts off at 1200 
keV. This peak corresponds to the radiation into the first harmonic of the ordinary 
channelling radiation (see e.g. p2|), and there is almost no radiation into higher 
harmonics. The latter fact is consistent with general theory of dipole radiation by 
ultra-relativistic particles undergoing quasi-periodic motion (see e.g. fll2f ). Dipole 
approximation is valid provided the corresponding undulator parameter p c = 27f-f(a c / A c ) 
is much less than 1. In this relation a c and A c stand for the characteristic scales of, 
correspondingly, the amplitude and the wave-length of the quasi-periodic trajectory. For 
the channelling motion one can estimate a c as d/2, and A c = cr c , where r c ~ 2ti J mj / U" 
standing for the characteristic period of the channelling oscillations (using the harmonic 
approximation for the interplanar potential one gets U" ~ 8Uo/d 2 where Uq is the 
depth of the interplanar potential well). In the case of 0.5 GeV positron channeled 
along (110) planes in Si one has 7 « 10 3 , U = 23 eV, d = 1.92 A [O. Hence, 



p c pa 0.2 <C 1 and all the channelling radiation is concentrated within some interval 
in the vicinity of the energy of the first harmonic. The latter one can estimate as 



(see e.g. [0) uj^ ~ Air^y 2 / V c ~ A^ 2 c/d^U /e arriving at the value huf^ pa 1190 keV 
which is exactly the calculated maximum value 1190 keV. More accurate estimates (the 
details are presented in |10||), based on the account for the dependence of the channeling 
oscillation period r c on the amplitude a c of the oscillations, also reproduce the calculated 
position and the width of the peak of the channeling radiation. 

Increasing the a/d ratio leads to the modifications in the spectrum of radiation. 
The changes which occur manifest themselves via three main features, (i) the lowering 
of the ordinary channelling radiation peak, (ii) the gradual increase of the intensity of 
undulator radiation due to the crystal bending, (iii) the appearing of additional structure 
(the sub-peaks) in the vicinity of the first harmonic of the ordinary channelling radiation. 
Let us discuss all these features. 

The decrease in the intensity of the ordinary channelling radiation with the increase 
of the a/d ratio is related to the simple fact that the growth of the AW amplitude leads 
to lowering of the allowed maximum value of the channelling oscillations amplitude a 
(this are measured with respect to the centerline of the bent channel) ]7J . Indeed, inside 
the bent channel the motion of the particle is subject to the action of the effective 



c 



9 



potential 

Ueff(p) = U(p)- 1 j^p (6) 

where p G [—d/2,d/2] is the (local) distance from the centerline and R(z) is the (local) 
curvature radius of the channel. For a channel bent by a transverse harmonic AW 
R = [(A/27r) 2 /a] sin(27rz/A). The particle could be trapped into the channelling mode 
provided its total energy associated with the transverse motion is less the minimal 
value, U e ff(p ), of the two maxima points of the asymmetric potential well described 
by @ 0, 0. The potential U e ff(p) reaches the magnitude of C/ e //(po) at some point 
p which satisfies the condition |p | < d/2 and the absolute value of p decreases with 
the growth of a. Hence, the larger the channel is bent the lower the allowed values 
of the channelling oscillations amplitude are, and, consequently, the less intensive is 
the channelling radiation, which is, essentially, proportional to the squared amplitude 



12fl . Let us note, that since we have restricted the range of the a c values by imposing 



the condition a c < d/2 — cltf (see the discussion above (|J)), then the decrease of the 
intensity of channelling radiation occurs starting with some non-zero value of the AW 
amplitude for which \p \ also becomes less than d/2 — a TF . 

These arguments explain the fact that the intensities of the ordinary channelling 
radiation for a/d = and a/d = 1 are much alike, while for a/d > 1 it starts loosing 
the magnitude. 

The undulator radiation (the AIR) related to the motion of the particle along the 
centerline of periodically bent channel is absent in the case of straight channel (the graph 
a/d = 0), and is almost invisible for comparatively small amplitudes of the AW (see the 
graph for a/d = 1). With a increasing the peaks corresponding to AIR are becoming 
more prominent, and for large a values (a/d ~ 10) two additional features appear: 
the intensity of the first harmonic of the AIR becomes larger than the intensity of the 
ordinary channelling radiation, and there appears radiation into the third harmonic of 
the AIR. 

The positions and the widths of the AIR peaks can be quite accurately estimated as 
follows (see [l], [|): = 87T7 2 cA _1 n/ (2+p 2 ), where the integer n = 1, 2 . . . enumerates 
the harmonics, and p = 27cy(a/X) stands for the parameter of the undulator related to 
the periodicity of the channel bending. The width of each peak Acu = (1/N)(uj^ /n) is 
independent on n. Substituting the values of 7, A and d indicated above one expresses the 
undulator parameter via the ratio a/d: p ph 0.05(a/d). Even for the largest considered 
value a/d = 10 the parameter p is less than 1, thus making the radiation into higher 
harmonics of the AIR almost negligible compared with the intensity radiated into the 
fundamental harmonic n = 1. The latter is located at %uj^ 1 ' ph 90 keV having the width 
hAu ph 6 keV which is almost 40 times less than the width of the peak of the channeling 
radiation. These values depend neither on the ratio a/d = 10 nor on the type of the 
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interplanar potential. 

As mentioned, all graphs in figure [3] refer to the forward emission. Therefore, in 
accordance with general theory of the undulator radiation (see e.g. |I2"||), the second 
peak of the AIR, which is mostly pronounced in the case a/d = 10, corresponds to the 
third harmonic of AIR, and is located at hto^ « 270 keV. The intensities radiated into 
the fundamental and the third harmonics are equal to 1.10 x 10 5 sr -1 and 7.12 x 10 3 sr _1 , 
respectively. Their ratio is approximately equal to p~ 4 which is also in accordance with 
general theory. 

The intensity of the AIR gradually increases with the AW amplitude. More precisely, 
di^AiR oc p 2 oc (a/d) 2 (in the case p < 1), and this tendency one can observe comparing 
the AIR peaks in the graphs corresponding to a/d = 4, 6, 8, 10. 

It is important to note that the positions of sharp AIR peaks, their narrow widths, 
and the radiated intensity are, practically, insensitive to the choice of the approximation 
used to describe the interplanar potential. In addition, provided the condition (0) is 
fulfilled, the AIR peaks are well separated (in the photon energy scale) from the peaks 
of the channeling radiation. Therefore, if being interested in the spectral distribution 
of the AIR only, one may disregard the channeling oscillations and to assume that the 
projectile moves along the centerline of the bent channel |l|, §. The above statements are 
illustrated by figure |] where we compare the results of calculations of the total spectrum 
(H) in vicinity of the first harmonic of AIR in the case a/d = 10. All parameters are 
the same as in figure |3|. The filled and open circles represent the results of evaluation of 
the right-hand sides of (Q) and ([5]) accompanied by numerical solution of the equations 
of motion for the projectile within the Moliere (filled circles) and the harmonic (open 
circles) approximations for the interplanar potential. The solid line corresponds to the 
AIR radiation only (in this case the numerical procedures are simplified considerably, 
leading to the reduction, by orders of magnitude, of the CPU time). It is clearly seen 
that the more sophisticated treatment almost does not change the profile of the peak 
obtained by means of simple formulae describing purely AIR radiation [jl], . Moreover, 
the minor changes in the position and the height of the peak can be easily accounted 
for within the framework of the cited formalism by introducing the effective undulator 
parameter |7|, 10 1 and (in the case of the harmonic approximation) the effective undulator 
amplitude [fT0|l . 

Thus, in vicinity of the AIR peaks there is no coupling of the two mechanisms of 
the radiation. On the contrary, the AIR strongly affects the spectrum in the photon 
energy range corresponding to the channeling radiation. The discussion of this effect 
lies beyond the scope of this Letter and is presented in [I0|. Here we only mention that, 
as it is clearly seen from the graphs a/d = 6, 8, 10 in figure ^, for large values of a/d the 
peaks of the ordinary channeling radiation acquire additional structure: there appear 
sub-peaks separated (in the case of the forward emission) by the interval 8oj = 2u^\ 
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Figure 4. Comparison of different approximations used to calculate the total 
ardiative spectrum in vicinity of the first harmonic of the AIR. The ratio a/d = 10, 
other parameters as in figure [| See also the commentaries in the text. 



The presented results of the numerical calculations of the total spectrum of radiation 
formed in an acoustically bent crystal clearly demonstrate the validity of the statement 
made in Jl], |2|, |F| that the AIR and the ordinary channeling radiation occur in essentially 
different ranges of the emitted photons energies, allowing, thus, to investigate the AIR 
properties separately from the ordinary channeling radiation. 

The positions of the peaks of the AIR as well as their intensity can be easily varied 
by changing the energy of projectile, the crystal type, the type of crystallographic plane, 
and the parameters of periodic pattern of crystal bendings. This statement by no means 
is restricted to the case of the acoustically bent crystal, but has more general nature. 
The treatment and the main results presented above can be applied, in particular, to 
consider the undulator radiation in statically bent crystals |2|, [|. Whatever way is 
chosen to prepare periodically bent crystalline lattice it is worth investigating, both 
theoretically and experimentally, not only the spontaneous emission of the AIR-type, 
but the stimulated emission as well [Q]. 
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